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AN ASCII FORTRAN SUBROUTINE FOR COMPUTING THE ELECTROMAGNETIC RADIATION 
SCATTERED BY A SPHERE 


INTRODUCTION 


The Lorenz-Mie formalism!s2:3 jis the analytical solution for the far 
field radiation which results when plane electromagnetic waves are scattered 
by a homogeneous spherical particle. With the advent of modern high speed 
computers programmed to provide numerical results, the Lorenz-Mie formalism 
may be regarded as an idealized experimental tool which is capable of 
predicting with certainty the scattering produced by spherical 
particulates. 


The development and use of computer codes to perform Lorenz-Mie 
calculations are of prime importance for smoke/aerosol research. For 
example, the excellent and widely used FORTRAN subroutine DBMIE,4 which 
was developed by J. V. Dave,> has been the mainstay for performing such 
calculations at Chemical Systems Laboratory (CSL). The DBMIE code was 
modified for use on the ARRADCOM UNIVAC 1108 computer as a FORTRAN V 
subroutine. The FORTRAN V compiler will be replaced in the near future 
by the UNIVAC 1100 series ASCII FORTRAN compiler.’ The ASCII compiler is 
a modern, state-of-the-art compiler which implements the ANSI FORTRAN 77 
standard. Rather than carry out a= straightforward conversion of the 
FORTRAN V Mie scattering subroutine, the code was completely rewritten to 
incorporate new ASCII FORTRAN language elements. The result was the ASCMIE 
subroutine which is a more structured, modular, Mie scattering code which 
executes faster than the FORTRAN V code. 


2 SOFTWARE DESIGN AND CODING CONVENTIONS 


A compilation listing of the ASCMIE subroutine is given in appendix A; 
the reader should become familiar with this listing before proceeding. This 
subroutine is an implementation of the downward recursion algorithm of 
Dave.o A call to subroutine ASCMIE requires the specification of the 
following input arguments: 


X= 7D/X, the size parameter for a sphere of diameter D 
illuminated by electromagnetic radiation of wavelength A. 
An error condition is generated and the execution of ASCMIE 
is halted if X < 0. X is a double precision variable. 


the real part of the complex refractive index, m. A 
printed warning is generated if N <1 because’ the 
algorithm has been thoroughly tested for N > 1 only. N is 
a double precision variable. 


imaginary part of the complex refractive index, m. An 
error condition is generated and the execution of ASCMIE is 
halted if K < 0. K is a double precision variable. 


THETAD = 6 double precision array containing the scattering angles 
In degrees. The maximum number of scattering angles Is set 
by the PARAMETER variable MXTHET at 100 angles. An error 
condition Is generated and the execution of ASCMIE is 
halted If any value of the scattering angle exceeds 90 
degrees. 
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NTHETA = the number of scattering angles at which the Lorenz-Mie 
solutions are sought. NTHETA should not exceed 100 unless 
ASCMIE Is recompiled with MXTHET set larger. An error 
condition ts generated and the execution of ASCMIE is 
halted If NTHETA > MXTHET. NTHETA Is an Integer variable. 


All the [Input arguments specified above are checked for error and warning 
conditions before a decision is made to proceed or halt execution. Each 
error or warning is described by a printed message. Also, a check is made 
to determine if sufficient array storage is available to perform the 
calculations before the program Is allowed to proceed. 


The subroutine returns the following outputs: 


QEXT = efficiency factor for total extinction. QEXT is a double 
precision variable. 
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QSCAT = efficiency factor for scattering. QSCAT is a double 
precision variable. 


, 


CTBRQS = product of the asymmetry factor and QSCAT. CTBRQS is a 
double precision variable. 


FMX = elements of the transformation matrix. FMX (4,MXTHET,2) is 
a double precision array. 


Definitions of the efficiency factors, asymmetry factor, and the 
transformation matrix can be found in any of several excellent treatises on 
the Lorenz-Mie Theroy.!> 3 


An attempt was made to make the ASCMIE code as structured as possible. Yq 
The beginning of this program module is indicated by the SUBROUTINE 
statement and the program module ends with the END statement. All of the 
FORTRAN code which implements the algorithm between the SUBROUTINE and ENO 
statements was divided into blocks according to the following scheme: 


(a) BLOCK 0 contains all specifications statements. ; 


(b) BLOCKS 100 to 600 are blocks of FORTRAN code which implement ‘| 
4 algorithmic tasks. 


(c) Each BLOCK is divided into as many SUBBLOCKS as necessary for 
clarity and convenience. 


ee 


Each BLOCK of the code was constructed according to the following style 
rules: 


(a) Each BLOCK starts with a statement labelled: 
0, 100, 200, ..., 600 
(b) Each SUBBLOCK starts with a statement labelled: 


O/106290 « ac lecye 290 
100,110,120, . . .,190 
200,210,220, . . .,290 


600,610,620, . . .,690 
The smallest available label is used. 


(c) 9000 series labels are FORMAT statements. Each digit following the 
9 serves to identify the BLOCK and SUBBLOCK in which the FORMAT 
statement occurs. 


Using the conventions described above, the ASCMIE code is represented by 
the following skeleton: 


START OF MODULE: SUBROUTINE ASCMIE (X, N, K, 
THETAD, NTHETA, QEXT, QSCAT, CTBROS, FMX) 


BLOCK 0: SPECIFICATIONS 

BLOCK 100: CHECK FOR INVALID OR 
OUT-OF-RANGE ARGUMENTS 

BLOCK 200: INITIALIZING CALCULATIONS 

BLOCK 300: SERIES SUMMATIONS 

BLOCK 400: CALCULATE FINAL RESULTS 

BLOCK 500: RETURN 

BLOCK 600: ERROR SECTION 

END OF MODULE: END 


3 =RESULTS 


3.1 Sample Problem. A FORTRAN main program, MIECHK, was written to call the 
ASCMIE subroutine and produce numeric results for sample problems. MIECHK 
is essentially the same program as that used previously by Dave? for the 
same purpose; a compilation listing of program MIECHK is given in appendix 
B. A MAP (Memory Allocation Processor) listing of the absolute element 
ASCMIE, which was produced by including MIECHK and ASCMIE, is given in 
appendix C. This absolute element was executed for Dave's sample problem 
m= 1.342 - 0.14, X = 1571.0.5 The results produced by ASCMIE are given 
in appendix D and are in agreement with previous results. 
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3.2 Timing Study. An important consideration in choosing a Mie scattering 
code is the computer time required to produce numeric results. The time 
required to execute ASCMIE for the sample problem described above (section 
3.1) was compared as the size parameter was varied from 0.1 to 1571 with 
similar results for two FORTRAN V versions of the Dave subroutine - MEMIE* 
and MIE2.** The results of the timing study are presented in table 1, where 
the times shown are the central processing unit times as listed by the 
UNIVAC 1108 accounting algorithm. 


Table 1. Comparison of Execution Times for Three Mie Scattering Subroutines 


sec 


0.1 
1.0 
10.0 
100.0 
1000.0 
1571.0 


In every case ASCMIE took less time to compute the same results thar the 
MEMIE and MIE2 subroutines. Independent of the value of X, we find that 


¥ aor (1) 
B where 
‘ Ta = time required when using the ASCMIE subroutine 
Ty, = time required when using the MEMIE subroutine 
Tug = time required when using the MIE2 subroutine 
Because of the uncertainties inherent in computer timing operations and 


the cursory nature of this timing study, these results should be regarded as 
suggestive but not definitive. 


*Unpublished data by M. E. Milham. 
**Unpublish -d data by . H. Frickel. 
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CONCLUSIONS 


The ASCMIE subroutine: 


(a) 
(b) 


(c) 
(d) 


Is based on the downward recursion scheme of Dave. 


Uses many of the FORTRAN 77 features available in the UNIVAC 1100 
series ASCII FORTRAN language. 


1s more modular and structured than the original Dave code. 


Executes faster than FORTRAN V versions of the Dave code (MEMIE, 
MIE2). 
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Research Division 
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DEPARTMENT OF DEFENSE 


Defense Technical Information Center 
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Director 
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ATTN: OB~4G1 

Washington, DC 20301 


Special Agent in Charge 
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Commander 
SED, HQ, INSCOM 


ATTN: IRFM~SED (Mr. Joubert) 
Fort Meade, MD 20755 


CEPARTMENT OF THE ARMY 


HQDA (DAMO-NCC) 
WASH DC 20310 


Deputy Chief of Staff for Research, 
Development & Acquisition 

ATTN: DAMA-CSS-C 

ATTN: DAMA-ARZ-0 

Washington, DC 20310 


US Army Research and Standardization 
Group (Europe) 

ATTN: DRXSN-E-SC 

Box 65, FPO New York 09510 


HQDA (DAMI-F1T) 
WASH, OC 20310 
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US Army Nuclear & Chemical Agency 
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Springfield, VA 22150 


Army Research Of fice 
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Research Triangle Park, NC 27709 
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US Army Medical Research and 
Development Command 
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Fort Detrick, MO 21701 
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USA Medical Research Institute of 
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NS Army 
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Project Manager Smoke/Obscurants 
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Aberdeen Proving Ground, MD 21005 
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US Army Foreign Science & Technology Center 
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US Army Missile Command 

Redstone Scientific Information Center 
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Redstone Arsenal, AL 35809 


Director 

DARCOM Field Safety Activity 

ATTN: ORXOS-C ! 
Charlestown, IN 47111 


Commander 


US Army Natick Research and 

Development Command 
ATTN: DRONA-VR \ 
ATTN: ORDNA-VT 1 


Natick, MA 01760 
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USA ARRADCOM 
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Director 


Ballistic Research Laboratory 
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ATTN: ORDAR-TSB-S 


Aberdeen Proving Ground, MD 21005 
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Readiness Command 

ATTN: ORSAR-ASN 

ATTN: ORSAR-POM 
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Commandant 
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Fort Benning, GA 31905 
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Commander 
USA Tralning and Doctrine Command 
ATTN: ATCD-Z 


Fort Monroe, VA 23651 


Commander 

USA Combined Arms Center and 
Fort Leavenworth 

ATTN: ATZL-CA-COG 

ATTN: ATZL-CAM-IM 

Fort Leavenwortt., KS 66027 


Commander 


US Army TRADOC System Analysis Activity 
ATTN: ATAA=SL 


White Sands Missile Range, NM 88002 


Commander/Director 
Atmospheric Sciences Laboratory 


ATTN: DELAS-EO (Dr. Niles) 
ATTN: DELAS-EO-MO 


White Sands Missile Range, NM 88002 
US ARMY TEST & EVALUATION COMMAND 


Commander 

US Army Test & Evaluation Command 
ATTN: DRSTE-CM-F 

ATTN: DORSTE-CT-T 

Aberdeen Proving Ground, MD 21005 


DEPARTMENT OF THE NAVY 


Commander 

Naval Air System Command 
AIR=310C (H. Rosenwasser) 
Washington, CC 20361 


Commander 
Naval Explosive Ordnance 
Disposa! Facility 
ATTN: Army Chemical Officer (Code AC-3) 
Indian Head, MD 20640 


Commander 
Naval Weapons Center 


ATTN: Technical Library (Code 343) 
China Lake, CA 93555 


Commanding Of ficer 

Naval Weapons Support Center 
ATTN: Code 5042 (Or. B.E. Douda) 
Crane, IN 47522 


US MARINE CORPS 


Clrector, Development Center 
Marine Corps Development and 
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Commander 
US Army Environmental Hygiene Agency 


Education Command ATTN: Librartan, Bldg 2100 
ATTN: Fire Power Division 1 Aberdeen Proving Ground, MD 21010 


Quantico, VA 22134 


Stimson Library (Documents) 1 
OEPARTMENT OF THE AIR FORCE Academy of Health Sciences 


HQ Foreign Technology Division (AFSC) 
ATTN: TQTR 1 


weight-Patterson AFB, OH 45433 


Fort Sam Houston, TX 78234 


HQ AFLC/LOWMM 1 
Wright-Patterson AFB, OH 45433 


OUTSIDE AGENCIES 


Battelle, Columbus Laboratories 
’ ATTN: TACTEC ! 


‘ 505 King Avenue 
| Columbus, OH 43201! 


Toxicology Information Center, WG 1008 
National Research Council 1 
2101 Constitution Ave., NW 

Washington, OC 20418 


